install.packages('aTSA')
library(aTSA)
data <- read.csv("data/时间序列分析——基于R（第2版）案例数据/csv/A1_7.csv")
ts_number <- ts(data$number)
# 时序图
plot(ts_number)
# 白噪声检验
# 纯随机性检验
for( k in 1:3) print(Box.test(ts_number, lag=6*k, type="Ljung-Box"))

# 平稳性检验
adf.test(ts_number)
